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ABSTRACT 

We investigate equilibrium sequences of magnetized rotating stars with four kinds of 
realistic equations of state (EOSs) of SLy (Douchin et al.), FPS (Pandharipande et al.), 
Shen (Shen et al.), and LS (Lattimer & Swesty). Employing the Tomimura-Eriguchi 
scheme to construct the equilibrium configurations, we study the basic physical prop- 
erties of the sequences in the framework of Newton gravity. In addition we newly take 
into account a general relativistic effect to the magnetized rotating configurations. 
With these computations, we find that the properties of the Newtonian magnetized 
stars, e.g., structure of magnetic field, highly depends on the EOSs. The toroidal mag- 
netic fields concentrate rather near the surface for Shen and LS EOSs than those for 
SLy and FPS EOSs. The poloidal fields are also affected by the toroidal configurations. 
Paying attention to the stiffness of the EOSs, we analyze this tendency in detail. In 
the general relativistic stars, we find that the difference due to the EOSs becomes 
small because all the employed EOSs become sufficiently stiff for the large maximum 
density, typically greater than lO^^g cm~^. The maximum baryon mass of the mag- 
netized stars with axis ratio g ~ 0.7 increases about up to twenty percents for that of 
spherical stars. We furthermore compute equilibrium sequences at finite temperature, 
which should serve as an initial condition for the hydrodynamic study of newly-born 
magnetars. Our results suggest that we may obtain information about the EOSs from 
the observation of the masses of magnetars. 
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1 INTRODUCTION 

Recently there is a growing evidence of the supermagnetized 
neutron stars with the mag netic fields of ~ lO" - 10^'^ G, 
the so-called magnet ars (e.g., IWoods fc ThompsonI (|2004h : 
iLattimer fc PrakashI (|2006l )). Although such stars are 
estimated to be only a subclass (~ 10% ) of the canonical 
neutr on stars with ~ lO'^^ — lO'^^ G (|Kouveliotou et al] 
1 19981 ). much attention has been paid to this objects because 
there remain many astrophysically exciting but unresolved 
problems. Giant flaring activities observed during the past 
~20 years have given us good opportunities to study the 
coupling of th e interior to the magnetospheric structures 
of magnetars (|Thompson fc DuncanI [l995l . Il996t ). but the 
relationship between the crustal fraction and the subs equent 
starqu arks has not yet been clarified (see references in lWattg 
(|2006l )). The origin of the large magnetic field is also a big 
problem, whether generated at post-collapse in the rapidly 
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neutron star (|Thompson fc DuncanI 19931) or descend from 



the main sequence stars ( Ferrario fc WickramasinghdbOOTl ) 



Assuming the large magnetic fields before core-collapse, 
extensive magnetohydrodynamic ( MHD) stellar collapse 



Kotake et al.l 2004 Moiseenko et al. 


20061: Takiwaki et all 


2004 lOberEauliuEer et al.l |2006|; 


Shibata et al.l l2006l: 


Livne et al.l 200?!) in order to understand the formation 



mechanism of magnetars with the different levels of the 
sophistication in the treatment of the equations of state 
(EO Ss), the neutr i no tra nsport, and the general relativity 
(see iKotake et al] (|2006l ) for a review). Here it is worth 
mentioning that the gravitational waves from magnetars 
could be a useful to ol giving u s an information about the 
magnetar interiors Jloka l200ll : iBonazzola fc GourgoulhonI 



1 19951 : iKotake et al] |2004| ). While in a microscopic point 
of view, the effects of magnetic field larger than the 
so-called QED limit o f ~ 4 .4 x 10^^ G, on the EOSs (e.g., 
ILattimer fc PrakashI (|2006l )) and the radi ation processes 
have been also extensively investigated (see iHarding fc Lail 
(|2006l ) for a review) . For the understanding of the formation 
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and evolution of the magnetars, it is indispensable to unify 
these macroscopic and microscopic studies together, albeit 
not an easy job. 

When one tries to study the above interesting issues, 
the construction of an equilibrium configuration of mag- 
netars may be one of the most fundamental problems, 
thus has be en extensively investigated sin c e the p i oneer - 
ing study byjchandrasekhar fc Fermil (|l953h . iFerrard (|l954l ') 
an d lRobertsllToSsI T studied the equilibrium con figurations of 
an inc ompressible star with the poloidal fields. iPrendergastI 
l|l956l ) then succeeded to take into account the toroidal 
fields. IWoltiej (|l960l) extended Prenderga st's w ork to the 
compr essible one. MonaghanI (|l965l . Il966l ) and iRoxburghl 
l| 19661 ) improved a boundary condition of the poloidal 
fields. Note that all these works mentioned above were 
based on the perturbation approach i n which one assumes 
the m agne tic field is sufficie ntly weak. lOstriker fc Hartwickl 
l| 19681 ) and lMiketinad (| 19751 ) developed techniques to treat a 
strong magnetic field be yond the perturbative a pproach with 



strong magnetic nem De yong tne pertur Dative a pproacn witn 
the N ewto nian gravity. Bocguet et al.l (|l995h , iKonno et al.l 
l|l999l ) and lloka fc Sasakil (|2003| . |2004| ) developed the scheme 
with the general relativistic framework. 

Here it should be noted that the studies mentioned 
above have c ontain e d two d r awbac k s to be modified. In 
iPrendergasj (Il956l). IWoltid (Il960l ). lOstriker fc Hartwickl 
( 19681 ) and lloka fc Sasakil (12004*) . the magnetic fields are 
assumed to vanish at the stellar surface and such a con- 
dition may be somewhat unrea l istic. I n iFerrarol 
RobertsI ([l95i ) . iMonaghanI (Il965l. Il966l). iMiketinad 



Bocguet et al 



(jl995l ) and lKonno et alf(| 19991 ) . onlv poloidal 



1954 ). 



19751 ). 



fields have been considered. It has been pointed out that 
such a purely poloidal field will be unstable and will decay 
within several Alfven times (IWrightlll973l:lMarkev fc Tavled 
1 19731 . 1 19741 ). In fact, iBraithwaite fc Spruit! (|2004l . |2006| ) 
recently find by the 3D numerical simulations that the 
pure poloidal fields certainly decay within several Alfven 
times and the configurations consisted both of poloidal and 
toroidal ones with comparable streng t h can be only stable. 

Recently, iTomimura fc Eriguchil (|2005l ) established a 
new non-perturbative method, which has conquered the 
above drawbacks. In the scheme, one can obtain a solution 
with a natural boundary condition, in which the magnetic 
field vanishes at the infinity. Moreover, one can treat both 
poloidal and toroidal field s simultaneously. More recently, 
lYoshida fc Eriguchil (|2006l ) calculated a number of sequences 
of magnetized rotating star and discussed their basic prop- 
erties. Extending the To mimura-Eriguc h i sch eme to treat 
the differential rotation, lYoshida et al.l (|2006l ) found that 
the magnetic fields can be in a twisted-torus equilibrium 
structure and discussed the universality of such structures. 

There may still rema in a room to be s ophis- 
ticated in their studie s ("Tomi mura fc Eriguc El l2005l : 
lYoshida fc Eriguchil l2006l : [Voshida et al.1 |2006| ) that the 
polytropic EOSs have been used for simplicity. In general, 
the central density of the neutro n stars is considered t o 
be higher than the nuclear density (|Bavm fc Pethic^l 19791 ). 
Since we still do not have an answer about the EOS in such 
a higher regime , many kinds of the n u clear EOSs have been 
proposed (e.g., iLattimer fc PrakashI (|2006h ) depending on 
the descriptions of the nuclear forces and the pos s ible a p- 
pearance of the exotic physics (e.g., iGlendennind (|200ll )). 
While the stiffness of the polytropic EOS has been treated 



to be globally constant inside the star, it should depend 
on the density in the realistic EOSs. Remembering that the 
equilibrium configurations are achieved by the subtle and lo- 
cal balance of the gravitational force, centrifugal force, the 
Lorentz force, and the pressure gradient, it is a nontrivial 
problem how the equilibrium configurations change for the 
realistic EOSs. 

Therefore the firs t purpose in this pap e r is t o ex- 
tend the studies of Tomimura fc Eriguchil (|2005h and 



lYoshida fc Eriguchil (|2006l ) to incorporate the realis- 
tic EOSs and discuss t he equilibrium properties . Four 
kinds of EOSs of SLy (iDouchin fc Haenself l200ll). EPS 
(iPandharipande fc Ravenhalll Il989l), Shen (IShen et al.1 



Il998al ). and Lattimer-Swestv (jLattimer fc Swestvl 1991 



adopted, which are often employed in the recent literatures 
of the MHD studies mentioned above. In contrast to the 
case of the polytropic EOS, a maximum density enters as 
a new parameter. At first, we set this value as comparable 
as the nuclear density of pnuc ~ 2.8 x lO^^g cm"'', because 
the general relativistic (GR) corrections are rather small 
at this density regime as will be mentioned. If a maximum 
density of star is much larger than the nuclear density, we 
must take account into a GR effect. However, the fully GR 
approach to the magnetized equilibrium configuration has 
not been establ i shed yet except for the purely poloidal fields 
(|Bocguet et al.l [19951 ). Therefore, we consider here a new 
approach to take into account a GR effect approximately 
and see their effects on the equilibrium configurations. This 
is the second purpose of this paper. Applying this method, 
we furthermore compute equilibrium sequences at finite 
temperature, which should serve as an initial condition 
for the hydrodynamic evolutionary studies of newly-born 
magnetars. 

This paper is organized as follows. In Section[2l we shall 
briefiy summarize the basic equations, the numerical scheme 
and the employed equations of state. Numerical results and 
their analysis are presented in Section [3| Summary and dis- 
cussion follow in Section |4| 



2 NUMERICAL SCHEME AND MODELS 
2.1 Basic equations 

The basic equations describing an equilibrium state of 
a perfect conductive fluid are summarize d according to 
iLovelace et al] (jl986h : lYoshida et al.1 (|2006l ): 



1. Maxwell's equations: 
VaB" = 47rp„ 

f. 4:7V 

aoeT—7 Ti -a 

e VfcSe = — J , 
c 

e^'^VtBe = 0, 
VaS" = 0, 



(1) 

(2) 

(3) 
(4) 



where we use Gaussian units, with the charge and current 
densities pe and both measured in electrostatic units. 
Here i?" , B", and c are the electric field, magnetic field, 
and the light velocity, respectively, Va stands for the 
covariant derivative and tate for the Levi-Civita tensor in 
fiat three-dimensional space. 
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2. Perfect conductivity equations: 

v'' 

K + eabe— B'=0, (5) 
c 

where v"^ denotes the fluid velocity. 

3. Mass conservation: 

Va(pi;") = 0, (6) 

where p is the mass density. 

4. Euler equations: 

pv''\/tVa + VaP + pVa$ - -€abej''B^ = 0, (7) 

C 

where p is the pressure and $ is the gravitational potential. 
Note that the last term on the left-hand side of equation 
(O represents the Lorentz force being exerted on the per- 
fectly conductive fluid, which can, due to equation ((2}, be 
rewritten in terms of _B" as 

-e^i^'B'' = ^{B'^tBa - B'VaBt), (8) 
c Atv 

5. Poisson equation for the gravitational potential: 

VVa* = AnGp, (9) 

where G is the gravitational constant. 

In order to obtain stationary states of magnetized stars, 
the following assumptions are adopted in this frame work:(l) 
The magnetic axis and the rotation axis are aligned. (2) The 
matter and the magnetic field distributions are axisymmet- 
ric about this axis. (3) The sources of the magnetic fields, 
i.e., the current distributions are confined to within the 
stars. For the sake of simplicity we require two additional 
conditions that are not necessary to have stationary solu- 
tions. (4) The meridional circulation of the fluid is neglected. 
Consequently, the fluid flow can be described by the angular 
velocity i} only, which is given by — i^ip"", where ip stands 
for the rotational Killing vector. Due to this assumptions 
and the axisymmetric assumption, the mass conservation 
equation Q is automatically satisfied. (5) The barotropic 
equation of state, p = p{p), is adopted. Although we employ 
realistic equations of state in this paper, this barotropic con- 
dition can be maintained as explained in subsection 12.21 

Under these assumptions, the basic equations for de- 
termining strictures of magnetized rotating stars can be de- 
rived. For axisymmetric configurations, the divergence-free 
magnetic fields are automatically satisfied by introducing ^I" 
of R and z as follows: 

(10) 



B' 



Rd: 



5^ = iA*. 
RdR 



(11) 



The function 4" is sometimes called the flux function. Here 
we have employed the cylindrical coordinates {R,cj>,z), in 
terms of which the line element and the Levi-Civita tensor 
are, respectively, given by 



as = QabClx ax 

(-R<t>z = ^• 



: dR^ + R^ dch^ + dz^ , 



Introducing the vector potential Aa, defined as 



(12) 
(13) 

(14) 



we can confirm that the flux function ^ is nothing but the 
component of the vector potential A^. Equation ([5} gives 
us the electric flelds as 



C 



(15) 



Note that — 0, because is an axisymmetric function. 
From equations ([3]) and (|15p the angular velocity must be 
expressed as 



n = n(l'). 



(16) 



This relation is so metimes called Ferraro's law of isorota- 
tion (|Ferrar oil 19371 ). 

The integrability conditions for equation (O result in 
the following relations: 



5<F 



i?2 



c 47r 



(17) 
(18) 



where 5* and /j. are arbitrary functions of ^P, and 

k(*) = 5''(*). (19) 

Here the prime denotes the differentiation with respect to 
^f. The functional relation 5'( ^) is sometimes called the dy- 
nami cal torque-free condition (|Chandrasekhadl 19561 : iMestell 

From equations ([2]) and (|18p we obtain 

_d_ f}_d^\ 9^ 
dRKRdRj^ dz^ 



-47rp7?^[i?^f7n' +^(*)] 

-K(v[>)S(vI/), (20) 



which is generalized Grad-Shafranov equation. The first in- 
tegral of the Euler equation can be expressed as 



/dP 1 /'* 

— = -R^n^ + J p{u)du + c, 



(21) 



where C is a co n stant of integration. Following 
iTomimura fc Eriguchil l|2005h . we convert partial dif- 
ferential equations ([9]) and (|20|l into integral equations, 
given by 



(22) 



, , , . , f KS + 47T{p+R'Qn')pR' . 7,3w„o^ 

A^{r)sm(f>= / sm (pd r {26) 

J 4:nR\r — f\ 

where all boundary conditions for the potentials $ and 'if, 
i.e., that the potentials are regular everywhere, are included 
in the integral expressions, and we need not consider them 
any further. Concerning the arbitrary function k(^'), since 
the current vector has to vanish outside the star as men- 
tioned before, the function k needs to vanish outside the 
star. 

Concerning the functions characterizing the magnetic 
field distributions, we choose the following forms: 



m(w) = 



•S'(") = , ° , (m - Umax)'''^^9{u - Umax), 
K + L 



(24) 
(25) 

k{u) = S' = a{u — UmaxtQ(u — U,nax) , (26) 

where p, a, and k are arbitrary constants and Umax is the 
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maximum value of *!/ in the vacuum region. Here 0{x) stands 
for the Heaviside's step function, and we have assumed that 
fc J5 in order to have finite k. It should be noted that, 
as shown below, k, defined by equation (|26p . satisfied the 
required condition, i.e., n — outside of the star, at least 
for the stars in the present study. 

Since we have not known what rotation law is realized 
in actual magnetized stars, in the present investigation we 
adopt the rigid rotation, 



r2(u) = Qo. 



(27) 



In order to carry out numerical computations properly, we 
introduce the following non-dimensional variables by using 
the maximum density pmax and the equatorial radius r^: 



p = pI pmax 



P "i-^Gplrax ' 
Cl = Q/ ^AnGpmax, 
C =C/{AnGpmaxr^e), 

- /' 1 

= ft/ 



M = M/ 



= A^/{V47TGpmaxr^), 
= Ba/{V4:nGpmaxre), 



(28) 

(29) 

(30) 
(31) 

(32) 

(33) 

(34) 
(35) 



where /3 is a numerical factor that is introduced to fix the 
non-dimensional equatorial radius to be unity. Here Ba de- 
notes the orthonormal component of Ba, which is frequently 
convenient because it has a physical dimension. By using 
these variables, global physical quantities characterizing an 
equilibrium state, the gravitational potential energy W , the 
rotational energy T, the internal energy U, the magnetic 
energy H, the mass M can be expressed as follows: 

W = W/(AnGpl,yt), (36) 

f = T/{4nGpl,y,), (37) 

U = U/{4nGpl,y,), (38) 

H = H/{4nGpl,,y:), (39) 

M = M/{pmaxrl), (40) 

where 
W ; 



2L 



-4 



u 



H = -- 



pimyd^r, 



pd^r, 



X eabej B d r, 



M ; 



/ 

J star 



pd^r. 



(41) 
(42) 
(43) 
(44) 
(45) 



the Hachisu self-consistent field ( HSCF 



1986 



We emplw 

scheme (iHachisul 
lYoshida fc Eriguchil 
HSCF scheme, one of the model parameters characterizing 
an equilibrium star is the axis ratio q, defined as g = rp/ve, 



iTomimura fc Eriguchil I2OO5I : 
20061 : lYoshida et al.1 I2OO6I I. In the 



where is the smallest distance to the surface from the 
origin. With the HSCF scheme, p, A^, f3, C, and f^o for 
rotating configurations (for non-rotating configurations) are 
iteratively solved, and during iteration cycles, q and other 
model parameters are fixed. By changing the axis ratio q, 
and fixing an appropriate set of the parameters, we follow 
one model sequence of equilibrium configurations. In actual 
calculations, we divide the interval [0, 1] in the r -direction 
into 100 meshes and the interval [0,7r/2]in the 9- direction 
into 200 meshes. Note that it is enough to calculate 
solutions for the interval [0,7r/2] in the ^-directions because 
we impose the equatorial plane symmetry. The accuracies 
of the numerical solutions are checked w ith an as s essme nt 
by the normalized virial equation (e.g., ICowlinj (|l965h ). 
defined as 



VG = \2T + W + 3U + H\/\W\. 



(46) 



For later convenience, here we like to explain about the 
qualitative meanings of the parameters a and fi. The param- 
eter fi is directly involved in the Bernoulli equation (|2H) and 
plays a role to determine the matter distribution. Increas- 
ing of a enhances magnitude of magnetic field ^ through 
the generalized Grad-Shafranov equation (|20|l . Therefore, 
increasing of both a and p. results in an enhancement in 
the Lorentz force exerting on the conductive fluids. Note 
that a — means that the magnetic field has the poloidal 
component only, however p, = does not mean that there 
exists only toroidal component (see equation (|10p . pi|l and 

m)- 



2.2 Equations of State 

As mentioned in Section [TJ equation of state (EOS) 
is an important ingredient for determining the equi- 
librium configurations. Before conducting an exten- 
sive study as done before for the studies of ro- 
tating equilibrium configurations in which mo r e var i- 
ety of EQSs w ere employed (e.g., iNozawa' et al.l (|l998h : 
iMorrison et al.l (|2004h and refe rences therein), we adop t 
here four kinds of EQSs of SLv (iDouchin fc Haense]||200ll) 



FPS jPandharipande fc Ravenhalllll989l). Shen jShen et al.1 
Il998al ). and Lattimer-Swestv (|Lattimer fc Swestvl Il99ll ) 
which are often employed in the recent MHD studies rel- 
evant for magnet ars. 

In the study of cold neutron stars, the /3-equilibrium 

condition with respect to beta decays of the form +p < > 

n + Us and n < > p + + Ue, can be well validated as for 

the static properties. Since neutrinos and antineutrinos es- 
cape from the star their chemical potentials vanish at zero- 
temperature T = with T being the temperature, and the 
equilibrium condition is pn ~ pc + Pp, with /in, pc, and pp 
being the chemical potentials of neutron, electron, and pro- 
ton, respectively. With the charge neutrality condition, we 
can determine the three independent thermodynamic vari- 
ables, (for example the pressure as P{p,Ye,T) with Ye be- 
ing the electron fraction), can only be determined by a single 
variable, which we take to be the density, namely P{p, Ye{p)) 
(|Shapiro fc Teukolskvl[l983l ). noting here that T = is as- 
sumed for the case of the cold neutron stars. Thanks to 
this, we can use the formalism mentioned in Section 2 with- 
out violating the barotropic condition of the EOSs. At the 
maximum densities higher than ~ 2/9nuc, muons can appear. 
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and h igher than ~ 3pnuc (|Wiringa et al.|[l988l : lAkmal et alj 
Il998|). one may tak e into possible appearance of hyperons 
( Glendennin3l200lh . However, since the muon contribution 
to pres sure at the higher density has been pointed to be very 
small (jPouchin fc Haense]||200ll ). and we still do not have 
detailed knowledge of the hyperon interactions, we prefer 
to employ the above neutron star matter, namely e~,n,p, 
model to higher densities. In the following, we shortly sum- 
marize features of the EQS s employed here. 

Lattimer-Swesty EOS (|Lattimer fc SwestvlflDQll ') is the 
one which has been used these years as a standard in 
the research field of core-collapse supernova explosions 
an d the subsequent neutr o n star formations (se e references 
in ISumivoshi et al.l (|2005l '): iKotake et al] Hooi)), which is 
based on the compressible drop model for nuclei together 
with dripped nucleons. The values of nuclear parameters 
are chosen according to nuclear mass formulae and other 
theoretical studies with the Skyrme interaction. Shen EOS 
IShen et all liggSal l is the rather modern one currently of- 
ten used in the research field, which is based on the rel- 
ativistic mean field theory with a local density approxi- 
mations, which has been constructed to reproduce the ex- 
perimental data of masses and radii of stable and unsta- 
ble nuclei (see references in IShen et al.l (|l998al )V FPS are 
m odern version of an earlier m i crosco pic EOS calculations 
bv [Friedman fc Pandharipandd (|l98lh . which employs both 
two body (U14) and three-body int eractions (TNI). In the 
SLy EOS (|Douchin fc Haensel|[200ll ) . neutron-excess depen- 
dence is added to the FPS EOS, which is more suitable for 
the neutron star interiors. As for the FPS and SLy EOSs, w e 
use the fitting formulae presented in IShibata et al. 
In the upper panel of Figure [TJ we plot the pressure P as 
a function of the rest-mass density p for the EOSs. In the 
bottom panel, we plot the adiabatic index, defined by 



dlnP 
dlnp ' 



(47) 



which can be a useful tool not only to characterize the stiff- 
ness but also to see their effects on the equilibrium configu- 
rations. As seen, the adiabatic index as a function of the den- 
sity is far away from constant as in the case of the polytropic 
EOS, showing an increasing (but sometimes zigzag) trend up 
to the several nuclear density (for example ~ 2pnuc of SLy) 
and decreases slowly at higher densities due to the inter- 
play of the density dependence of the nuclear interactions 
and of increasing proton fractions. Note that the relatively 
large discrepancies of the pressure below ~ 10^* g cm""^ 
and the zigzag features of the adiabatic indices are due to 
the difference in the treatment of the inhomogeneous matter 
consisted of electrons, protons, free nucleons, and a species 
of heavy nuclei, and these are pointed out to have little ef- 
fects on determining the equilibrium configurations which 
predominantly d etermined by the behavior of the EOSs at 
higher densities (|Shen et al.lll998bh . 
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Figure 1. Pressure (top) and the adia batic index F (bo ttom) as a 
functio n of rest-mass density for SLy llDouchin fc HacIlscl^ 200^). 
FPS l lPandharipande fc Ravenhalll 19891') . Shen l lShen et all 
Il998al) . and LS llLattimer fc Swestvl Il99lj) (LS) EOSs, respec- 



tively. The vertical blue (short-dashed) and green (dashed) lines 
indicate the densities of 3 X lO'^* g cm~^ (Section 13.11 1 and 
1 X 10^^ g cm~^ (Section I3.2|l set hero for the maximum den- 
sities. 



this choice makes the maxim um toroidal fields to be c om- 
parable to the poloidal ones (|Yoshida fc EriguchillioO^ . It 
should be mentioned that one should specify the maximum 
density when using the realistic EOSs, while the degree of 
freedom of this choice can be eliminated in case of the poly- 
tropic EOSs because one can choose the polytropic constant 
freely. We actually fix its value as 3 x 10^* g cm^"^ close to 
the nuclear density. 

Before going to the main results, we note that we have 
checked the accuracy of our newly developed code by the 
test problems and their results are summarized in appendix 



in 



3 NUMERICAL RESULTS 

As mentioned, the first purpose of this study is to investi- 
gate the basic properties of magnetized rotating stars with 
the four kinds of EOSs in Newtonian gravity. In all the calcu- 
lations, the value of k (equation (|25p ') is set to be 0.1 because 



3.1 Newtonian Case 

For clarity, we categorize the equilibrium configurations to 
two types of non-rotating type and rotating one. For non- 
rotating sequence, the parameters, axis ratio q and a, are 
input and kept fixed during calculation and then the pa- 
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Figure 2. H/\W\ as function of q for the non rotating sequences 
shown in Table [T] 

rameters, /5, C, and fi are obtained as output. For rotating 
sequence, q, a, and fi are input and kept fixed, then /3, C, 
and f^o are obtained as outcome. 



3.1.1 Static Magnetized Configurations 

We first concentrate on the non-rotating sequences, in which 
the anisotropic magnetic stress is only the agent to deform 
the stars. In Table[T] various p hysical quantities correspond- 
ing to the previous studies (|Tomimura fc Eriguchil l2005l : 
lYoshida fc EriguchilbOoH : lYoshida et al.ll2006l ) are given for 
the four kinds of EOSs. As a model parameter, we set 
a = 12 because we are interested in the combination effects 
of poloidal and toroidal fields. Note again that the choice of 
a = leads to the unstable pure poloidal configurations as 
mentioned. 

As shown in the table, the values of _ff/|H^| for the se- 
quences become large enough to be 0.1 when the stars suf- 
ficiently deform typically q of the axis ratio to be < 0.7. In 
such a region, the perturbative approach may break down 
and a non-perturbative taken here is valid. 

Lowering the axis ratio, it can be seen that all sequences 
can achieve nearly toroidal density configuration, q ~ due 
to the strong Lorentz forces. In fact, iJ/j W| increases as the 
axis ratio q decreases as shown in Figure[21 This feature does 
not depend on the EOSs. Furthermore we find the values of 
for SLy or FPS sequences is greater than the those 
for Shen or LS sequences for any axis ratio. This tendency 
is explained with respect to the stiffness of the EOSs. As 
seen from the bottom panel of Figure 1, Shen or LS EOSs 
are more stiffer than SLy or FPS EOSs in the density region 
near and below 3x 10^* g cm~'^. Therefore, if the stars should 
have the same axis ratio, the force driven by the pressure 
for the Shen or LS stars is greater than the forces for the 
SLy or FPS stars. Consequently, the SLy or FPS stars need 
much Lorentz force than the Shen or LS stars to have the 
same degree of the deformation. 



3.1.2 Rotating Magnetized Configurations 

Tables [2] and [3] are devoted to the magnetized sequences 
with rotation for the four kinds of EOSs. By setting a = IQ 
and jl = 0.2 in Table [1 a = 16 and = 0.3 in Table O 



respectively, we can compute the equilibrium configurations 
with comparable toroidal and poloidal fields. By changing 
jl in two values, we can see clearly the effects of stronger 
magnetic field. 

From the tables, it can be seen that there exist maxi- 
mum {qmax) and minimum (gmin) values in the axis ratio q. 
This feature does not appear in the sequence only with the 
magnetic fields (see Subsection 13 . 1 . 1 1 ) . in which we can find 
the physical solutions for any values of q (see Table [l|. 

Equilibrium configurations which have larger axis ra- 
tio than qmax have < and such configurations are of 
course unphysical. Too strong magnetic field causes such a 
configuration, in which "anti" -centrifugal force is needed. On 
the other hand, configurations with smaller q than qmin be- 
long to following two types: (1) Due to too strong Lorentz 
forces, the converged solutions have Q? < as explained 
above. (2) the mass of the star sheds from the equatorial 
surface because the centrifugal forces are too strong (MS in 
the table indicates the so-cal le d mass-shedding). F ol lowin g 
Tomimura fc Eriguchil (|2005l ): lYoshida fc Eriguchil (|2006h : 



Yoshida et al.1 (|2006l ). we call the sequence ending in the for- 



mer and the latter type as magnetic-field-dominated (MD) 
and rotation-dominated (RD) sequence, respectively. It is 
noted that larger values of jl make qmax smaller (compare 
Table [5] with 13)) because qmax is determined by the magnetic 
field strength. 

So, all the sequences shown in Table [2] belong to the RD 
sequence regardless of the EOSs, but not for the sequences 
in Tabled The sequences with SLy or FPS EOSs belong to 
RD type and those with Shen or LS EOSs belong to MD 
type. As like this, even if we set the same parameter of a 
and jl, it is found that the type of the sequence depends on 
the EOSs. 

In order to see this feature more clearly, we perform a 
parameter search m a — ji parameter space and classify the 
sequences. In Figure (3] we show a set of phase diagrams on 
a — jl plane with different EOSs. In the figure, the region 
with < shows that we cannot obtain any solutions 
except one with negative angular velocity, "qmin = 0.01" 
indicates a parameter region in which nearly toroidal con- 
figurations exist. So we find the equilibrium sequences with 
SLy or FPS EOSs are classified into either RD or MD type. 
However, the sequences with Shen or LS has another type of 
sequence, in which nearly toroidal configurations g ~ exits 
as in the case of magnetized configurations without rotation 
(see Section l3.1.ip . Such a configuration never appears in 
the model with SLy or FPS EOSs. Looking at the stiffness 
of EOSs in Figure [T] again, it can be seen that the Shen and 
LS EOSs are stiffer than SLy and FPS EOSs near the cen- 
tral density adopted here [pmax = 3 x lO^^g cm~^). Thus 
the sequence with a nearly toroidal configuration is found 
to appear for the stiffer EOSs. This qualitative feature was 
also noticed in the polytropic studies (flbmimura fc Eriguchil 



also noticed m tne polytropic studies (I lomimura &i H 
I2OO5I : lYoshida fc Eriguchilbooel : lYoshida et al.ll200^ 



Then we move on to investigate the structures of the 
equilibrium configurations. Figures 1 4171 show the distribu- 
tions of density, and toroidal/poloidal magnetic fields of 
models characterized by g = 0.7, a — 16, ji = 0.2 for SLy, 
FPS, Shen, and LS EOSs, respectively. The reason why these 
values of the parameters are chosen is that we want to fix q 
for the different EOS models. For example, there does not 
exist the common value of q if we set the value of fi as 0.3 
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Table 1. Physical quantities for the sequence with a = 12 and = 0. 



q 


H/\W\ 


u/\w\ 




\w\ 






A 




C 






,8 




M 




VC 




SLy 


0.99 


0.004 


0.332 


0, 


140E- 


■01 


0, 


,079 


-0, 


.311E- 


■01 


0, 


,256E- 


■01 


0.387 


0, 


.125E- 


■03 


0.90 


0.040 


0.320 


0, 


168E- 


■01 


0. 


,224 


-0, 


,374E- 


■01 


0, 


,259E- 


■01 


0.431 


0, 


,130E- 


■03 


0.80 


0.091 


0.303 


0, 


,212E- 


■01 


0. 


,288 


-0, 


,468E- 


■01 


0, 


,260E- 


■01 


0.495 


0, 


,112E- 


■03 


0.70 


0.151 


0.283 


0, 


,277E- 


■01 


0, 


,319 


-0, 


.597E- 


■01 


0, 


,256E- 


■01 


0.582 


0, 


.828E- 


■04 


0.60 


0.212 


0.263 


0, 


,376E- 


■01 


0, 


,334 


-0, 


.772E- 


■01 


0, 


,243E- 


■01 


0.700 


0, 


.690E- 


■04 


0.50 


0.268 


0.244 


0, 


,460E- 


■01 


0, 


,342 


-0, 


.935E- 


■01 


0, 


,218E- 


■01 


0.795 


0, 


.813E- 


■04 


0.40 


0.328 


0.224 


0, 


,377E- 


■01 


0. 


,342 


-0, 


,905E- 


■01 


0, 


,184E- 


■01 


0.733 


0, 


,755E- 


■04 


0.30 


0.353 


0.216 


0, 


,245E- 


■01 


0, 


,309 


-0, 


.739E- 


■01 


0, 


,148E- 


■01 


0.600 


0, 


.105E- 


■03 


0.20 


0.359 


0.214 


0, 


182E- 


■01 


0, 


,290 


-0, 


.639E- 


■01 


0, 


,127E- 


■01 


0.520 


0, 


.149E- 


■03 


0.10 


0.364 


0.212 


0. 


154E- 


■01 


0, 


,282 


-0, 


,591E- 


■01 


0, 


,116E- 


■01 


0.479 


0, 


,143E- 


■03 


0.01 


0.366 


0.211 


0. 


145E- 


■01 


0. 


,280 


-0, 


,577E- 


■01 


0, 


,112E- 


■01 


0.466 


0, 


,133E- 


■03 



EPS 



0, 


,99 


0, 


,004 


0, 


,332 


0, 


,394E- 


■02 


0, 


,104 


-0, 


.145E- 


■01 


0, 


,155E- 


■01 


0, 


,180 


0, 


.231E- 


■03 


0, 


,90 


0, 


,040 


0, 


,320 


0, 


,577E- 


■02 


0, 


,282 


-0, 


.198E- 


■01 


0, 


,170E- 


■01 


0, 


,226 


0, 


.199E- 


■03 


0, 


,80 


0. 


,095 


0, 


,302 


0, 


,916E- 


■02 


0. 


,344 


-0, 


,286E- 


■01 


0, 


,187E- 


■01 


0, 


,299 


0, 


,161E- 


■03 


0. 


,70 


0. 


,166 


0, 


,278 


0, 


152E- 


■01 


0, 


,357 


-0, 


,424E- 


■01 


0, 


,199E- 


■01 


0, 


,406 


0, 


,114E- 


■03 


0, 


,60 


0. 


,238 


0, 


,254 


0, 


,256E- 


■01 


0, 


,351 


-0, 


.623E- 


■01 


0, 


,202E- 


■01 


0, 


,557 


0, 


.748E- 


■04 


0, 


,50 


0, 


,298 


0, 


,234 


0, 


,353E- 


■01 


0, 


,348 


-0, 


.809E- 


■01 


0, 


,187E- 


■01 


0, 


,682 


0, 


.632E- 


■04 


0, 


,40 


0. 


,361 


0, 


,213 


0, 


,277E- 


■01 


0. 


,338 


-0, 


,763E- 


■01 


0, 


,156E- 


■01 


0, 


,618 


0, 


,678E- 


■04 


0, 


,30 


0. 


,373 


0, 


,209 


0, 


163E- 


■01 


0. 


,292 


-0, 


,585E- 


■01 


0, 


,120E- 


■01 


0, 


,482 


0, 


,105E- 


■03 


0, 


,20 


0, 


,381 


0, 


,206 


0, 


121E- 


■01 


0, 


,277 


-0, 


.509E- 


■01 


0, 


,102E- 


■01 


0, 


,418 


0, 


.129E- 


■03 


0, 


,10 


0, 


,386 


0, 


,205 


0, 


102E- 


■01 


0, 


,271 


-0, 


.471E- 


■01 


0, 


,931E- 


■02 


0, 


,385 


0, 


.139E- 


■03 


0. 


.01 


0. 


,388 


0. 


,204 


0, 


,966E- 


■02 


0. 


,270 


-0, 


,459E- 


■01 


0, 


,903E- 


■02 


0, 


,374 


0, 


,145E- 


■03 



Shen 



0.99 


0.002 


0.332 


0.161 


0.051 


-0, 


132 


0.701E-01 


1.65 


0.218E-02 


0.90 


0.027 


0.324 


0.149 


0.157 


-0, 


134 


0.645E-01 


1.58 


0.207E-02 


0.80 


0.059 


0.313 


0.137 


0.220 


-0, 


137 


0.580E-01 


1.50 


0.197E-02 


0.70 


0.096 


0.301 


0.127 


0.266 


-0, 


141 


0.510E-01 


1.43 


0.185E-02 


0.60 


0.138 


0.287 


0.120 


0.301 


-0, 


147 


0.436E-01 


1.38 


0.169E-02 


0.50 


0.185 


0.271 


0.111 


0.328 


-0, 


151 


0.365E-01 


1.32 


0.156E-02 


0.40 


0.234 


0.255 


0.091 


0.345 


-0. 


146 


0.308E-01 


1.20 


0.140E-02 


0.30 


0.279 


0.240 


0.070 


0.348 


-0, 


134 


0.262E-01 


1.05 


0.135E-02 


0.20 


0.303 


0.232 


0.054 


0.338 


-0, 


120 


0.230E-01 


0.932 


0.144E-02 


0.10 


0.305 


0.231 


0.048 


0.325 


-0, 


113 


0.215E-01 


0.878 


0.146E-02 


0.01 


0.306 


0.231 


0.046 


0.322 


-0, 


110 


0.210E-01 


0.860 


0.148E-02 



LS 



0.99 


0.003 


0.331 


O.lOOE+00 


0.054 


-0, 


,010 


0, 


.570E- 


■01 


1.25 


0, 


,357E-02 


0.90 


0.031 


0.322 


0.965E-01 


0.166 


-0. 


104 


0, 


.532E- 


■01 


1.22 


0, 


,313E-02 


0.80 


0.067 


0.310 


0.931E-01 


0.229 


-0. 


110 


0, 


.485E- 


■01 


1.19 


0, 


,295E-02 


0.70 


0.108 


0.296 


0.910E-01 


0.273 


-0, 


117 


0, 


.434E- 


■01 


1.17 


0, 


,277E-02 


0.60 


0.154 


0.281 


0.909E-01 


0.307 


-0, 


126 


0, 


.377E- 


■01 


1.18 


0, 


,221E-02 


0.50 


0.203 


0.265 


0.887E-01 


0.331 


-0, 


134 


0, 


.320E- 


■01 


1.16 


0, 


,209E-02 


0.40 


0.255 


0.248 


0.738E-01 


0.345 


-0, 


130 


0, 


.270E- 


■01 


1.07 


0, 


,184E-02 


0.30 


0.299 


0.233 


0.553E-01 


0.343 


-0, 


118 


0, 


.229E- 


■01 


0.928 


0, 


,221E-02 


0.20 


0.316 


0.227 


0.428E-01 


0.327 


-0, 


105 


0, 


.201E- 


■01 


0.822 


0, 


,202E-02 


0.10 


0.317 


0.227 


0.372E-01 


0.314 


-0. 


,098 


0, 


.187E- 


■01 


0.769 


0, 


,213E-02 


0.01 


0.318 


0.227 


0.355E-01 


0.311 


-0, 


,096 


0, 


.182E- 


■01 


0.751 


0, 


,191E-02 



(see Table ISj. q = 0.7 taken here is the smallest common 
value in the sequences with jl — 0.2 (see Table [2}. 

It is seen that the density distributions of stars with SLy 
or FPS EOSs are more concentrated at the center than those 
of the stars with Shen or LS EOSs. These feature are due to 
the stiffness of EOSs as mentioned above. We also find the 
toroidal fields for SLy or FPS EOSs distribute in relatively 
wider regions in the vicinity of the equatorial plane than for 



Shen or LS EOSs, in which the toroidal fields concentrate 
near the stellar surfaces. It is noted that the toroidal fields 
only exist in the interior of the stars, which is due to the 
choice of the functional form of k (see equation (|26|) ). Near 
the rotational axis, the poloidal magnetic field behaves like 
a dipole field, where the toroidal fields are weak. However 
in the vicinity of the equatorial plane where the toroidal 
fields become strong, the poloidal fields become distorted. 
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Table 2. Physical quantities for rotating magnetized stars with a = 16 and fi = 0.2. 



q 


H/\W\ 


U/\W\ 


T/\W\ 


\W\ 




C 


/3 


M 


h 


VC 


SLy 


Q2 < 






















0.90 


0.350E-1 


0.321 


0.855E-3 


0.147E-1 


0.536E-3 


-0.347E-1 


0.249E-1 


0.399 


0.320 


0.128E-3 


0.80 


0.283E-1 


0.316 


0.117E-1 


0.898E-2 


0.720E-2 


-0.290E-1 


0.202E-1 


0.297 


0.270 


0.167E-3 


0.70 


0.198E-1 


0.315 


0.176E-1 


0.452E-2 


0.108E-1 


-0.220E-1 


0.154E-1 


0.197 


0.210 


0.216E-3 


0.63 


0.146E-1 


0.316 


0.184E-1 


0.251E-2 


0.114E-1 


-0.172E-1 


0.122E-1 


0.138 


0.199 


0.261E-3 


MS 






















FPS 


^2 < 






















0.94 


0.168E-1 


0.328 


0.374E-3 


0.417E-2 


0.247E-3 


-0.157E-1 


0.156E-1 


0.187 


0.189 


0.225E-3 


0.90 


0.152E-1 


0.327 


0.234E-2 


0.341E-2 


0.154E-2 


-0.146E-1 


0.143E-1 


0.166 


0.180 


0.242E-3 


0.80 


0.113E-1 


0.326 


0.583E-2 


0.191E-2 


0.383E-2 


-0.116E-1 


0.114E-1 


0.117 


0.166 


0.306E-3 


0.70 


0.783E-2 


0.326 


0.729E-2 


0.941E-3 


0.480E-2 


-0.866E-2 


0.858E-2 


0.765E-1 


0.160 


0.402E-3 


0.65 


0.644E-2 


0.326 


0.742E-2 


0.632E-3 


0.489E-2 


-0.735E-2 


0.733E-2 


0.603E-1 


0.168 


0.477E-3 


MS 






















Shcn 


Q2 < 






















0.80 


0.605E-1 


0.311 


0.216E-2 


0.132 


0.138E-2 


-0.134 


0.573E-1 


0.146E+1 


0.424 


0.197E-2 


0.70 


0.672E-1 


0.289 


0.315E-1 


0.106 


0.183E-1 


-0.129 


0.485E-1 


0.129E+1 


0.458 


0.182E-2 


0.60 


0.756E-1 


0.265 


0.644E-1 


0.814E-1 


0.335E-1 


-0.122 


0.397E-1 


O.lllE+1 


0.497 


0.166E-2 


0.50 


0.860E-1 


0.237 


0.101 


0.596E-1 


0.461E-1 


-0.113 


0.309E-1 


0.929 


0.537 


0.148E-2 


0.40 


0.991E-1 


0.204 


0.144 


0.406E-1 


0.550E-1 


-0.101 


0.224E-1 


0.750 


0.582 


0.124E-2 


MS 






















LS 


Q2 < 






















0.81 


0.669E-1 


0.309 


0.143E-2 


0.881E-1 


0.883E-3 


-0.106 


0.482E-1 


0.115E-f-l 


0.502 


0.289E-2 


0.80 


0.675E-1 


0.307 


0.400E-2 


0.862E-1 


0.246E-2 


-0.106 


0.475E-1 


0.114E+1 


0.505 


0.291E-2 


0.70 


0.749E-1 


0.287 


0.313E-1 


0.675E-1 


0.176E-1 


-0.101 


0.401E-1 


0.988 


0.535 


0.275E-2 


0.60 


0.834E-1 


0.264 


0.615E-1 


0.495E-1 


0.312E-1 


-0.934E-1 


0.326E-1 


0.826 


0.565 


0.242E-2 


0.50 


0.905E-1 


0.240 


0.929E-1 


0.304E-1 


0.422E-1 


-0.799E-1 


0.246E-1 


0.622 


0.584 


0.238E-2 


0.48 


0.865E-1 


0.239 


0.972E-1 


0.245E-1 


0.439E-1 


-0.732E-1 


0.226E-1 


0.547 


0.567 


0.236E-2 


MS 























The region where the magnetic field is mixed depends on 
the EOSs. The stars with the softer EOS (SLy or FPS) tend 
to have wider mixed region than those with the stifi'er EOS 
(Shen or LS). The magnetic field lines 'i' may be good a tool 
to see the structures, showing the tori of twisted field lines 
around the symmetry axis inside the star and the untwisted 
poloidal field lines, which penetrate the stellar surface to 
continue to the vacuum. It is noted that this universality of 
the twisted-torus structures of the magnetic fiel ds was seen 
also in the polytropic stars (jYoshida et al.ll2006h . Therefore 
our results for the realistic EOSs may be regarded as a fur- 
ther generalization of their results. 

3.2 General Relativistic Case 

As mentioned in Section [1] when one tries to perform a 
qualitative investigation on the equilibrium configurations 
of magnetized rotating stars, the general relativistic effect 
must be taken into account. Since a full treatment is be- 
yond scope of our paper, we here try to include the effect 
by adding a spherically general relativistic correction to the 
Newtonian gravitational potential. Although this method is 



similar to th e appr oach reported in lRampp fc Jarikal (|2002l ): 
iMarek et al.l (|2006l ). the definition of general relativistic po- 
tential and its boundary condition are modified to be appro- 
priate for this study. Details of our approach and its verifica- 
tion are shown in appendix[B] In this subsection, we report 
the models with the correction. 



3.2.1 Basic property 

In Table U physical quantities of magnetized stars, such as 
magnetic field at pole Bp, one at center Be, baryon mass 
M, rotation period P, and radius R, are shown with a = 20 
and fi = 0.1. Maximum density is taken to be lO^^g cm^'^, 
typically where the general relativistic correction cannot be 
negligible. Since our treatment for the correction may break 
down if the equilibrium star is highly deformed, we only 
focus on the equilibrium sequences with mildly strong mag- 
netic field with the comparable strength of poloidal and 
toroidal fields (see h in Table UJ. 

From the table, it is found that there exist Qmax and 
Qmin as same in the Newtonian sequences. We find that all 
the sequences belong to the rotation-dominated type, irre- 
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Table 3. Same as Tabic [2] but with a = 16 and fl = 0.3. 



q 


H/\W\ 


u/\w\ 


T/\W\ 


\W\ 




C 


/3 


M 


h 


VC 


SLy 


Q2 < 






















0.70 


0.147 


0.282 


0.283E-2 


0.195E-1 


0.144E-2 


-0.494E-1 


0.232E-1 


0.474 


0.471 


0.820E-4 


0.60 


0.425E-1 


0.305 


0.214E-1 


0.323E-2 


0.125E-1 


-0.206E-1 


0.128E-1 


0.161 


0.203 


0.241E-3 


MS 






















FPS 


Q2 < 






















0.86 


0.492E-1 


0.317 


0.301E-3 


0.525E-2 


0.187E-3 


-0.195E-1 


0.162E-1 


0.214 


0.217 


0.207E-3 


0.80 


0.372E-1 


0.318 


0.483E-2 


0.326E-2 


0.304E-2 


-0.158E-1 


0.135E-1 


0.161 


0.179 


0.255E-3 


0.70 


0.226E-1 


0.321 


0.793E-2 


0.135E-2 


0.507E-2 


-0.107E-1 


0.967E-2 


0.950E-1 


0.159 


0.361E-3 


0.63 


0.161E-1 


0.323 


0.819E-2 


0.700E-3 


0.529E-2 


-0.802E-2 


0.750E-2 


0.640E-1 


0.167 


0.459E-3 


MS 






















Shcn 
























0.52 


0.190 


0.269 


0.518E-3 


0.106 


0.254E-3 


-0.146 


0.367E-1 


0.129E+01 


0.491 


0.153E-2 


0.50 


0.195 


0.264 


0.635E-2 


0.103 


0.304E-2 


-0.146 


0.350E-1 


0.127E+01 


0.496 


0.154E-2 


0.40 


0.222 


0.237 


0.332E-1 


0.850E-1 


0.136E-1 


-0.142 


0.275E-1 


0.115E-I-01 


0.525 


0.135E-2 


0.30 


0.260 


0.216 


0.459E-1 


0.654E-1 


0.151E-1 


-0.131 


0.222E-1 


O.lOlE-l-01 


0.565 


0.124E-2 


0.20 


0.305 


0.213 


0.273E-1 


0.505E-1 


0.713E-2 


-0.117 


0.201E-1 


0.897 


0.610 


0.139E-2 


0.11 


0.331 


0.222 


O.llOE-2 


0.457E-1 


0.259E-3 


-0.110 


0.200E-1 


0.858 


0.615 


0.148E-2 


< 






















LS 


Q2 < 






















0.51 


0.217 


0.259 


0.243E-2 


0.813E-1 


0.112E-2 


-0.128 


0.311E-1 


O.lllE-f-Ol 


0.569 


0.205E-2 


0.50 


0.220 


0.256 


0.512E-2 


0.805E-1 


0.233E-2 


-0.128 


0.304E-1 


O.llOE-fOl 


0.571 


0.194E-2 


0.40 


0.251 


0.230 


0.289E-1 


0.695E-1 


0.113E-1 


-0.127 


0.241E-1 


0.102E-I-01 


0.600 


0.171E-2 


0.30 


0.298 


0.212 


0.322E-1 


0.533E-1 


0.981E-2 


-0.117 


0.199E-1 


0.905 


0.648 


0.191E-2 


0.21 


0.350 


0.216 


0.842E-3 


0.418E-1 


0.204E-3 


-0.104 


0.186E-1 


0.811 


0.690 


0.105E-2 


< 























spective of EOSs, with the maximum and minimum value 
of a:xis ratio being almost the same, typically qmax ~ 0.9 
and Qmin ~ 0.6. More interestingly, we find the GR effect 
does not change the basic property of the magnetized stars 
drastically. 

The values of the magnetic field at pole and center are 
also comparable, typically Bp ~ 10^^ [G] and Be ~ 10"[G]. 
The baryon masses of these stars are in the range between 
1.7M0 - 3.6M0 and the stars with FPS EOS tend to have 
lighter masses and those with Shen EOS. Masses of the stars 
with SLy and LS EOSs are shown to be c omparable. C om- 
paring with the observations of magnetars (|Kaspill2004 ). the 
obtained magnetic fields are a bit larger and the rotation 
periods are an order of milliseconds, which are too rapid in 
comparison with the observed ones of seconds. We discuss 
this point in the Subsection 13.31 

Figures [S] HI 1101 and [TT] are the density and magnetic 
field distributions of the equilibrium configurations with 
q — 0.7 in Table |4l Note again that this value is the small- 
est common value in the configurations with the different 
EOSs. It can be seen that the configurations of density and 
the magnetic fields are rather similar regardless of the EOSs. 
The toroidal components of the magnetic fields concentrate 
near the stellar surface. The poloidal fields are distorted in 
the region, in which the toroidal fields are strong, and its 



shapes are dipole like near the rotation axis. Despite of the 
incursion of the general relativistic correction, the structures 
of the magnetic field lines are twisted tori as same in the 
Newtonian case. The reason of these similarity is that the 
employed EOSs become sufficiently stiff (Figure [1} below 
the maximum density of lO^^g cm~'' and that the effects 
on the configurations become smaller. It is seen that Shen 
and LS EOSs become softer at the density regime due to the 
increasing proton fraction. However it is noted that the equi- 
librium configurations are affected whether or not the EOSs 
become already stiffened at the smaller density regimes than 
the maximum density. In the two EOSs, this indeed occurs 
at ~ lO^^g cm~"^. After the stiffening, the effects of the dif- 
ference of the EOSs on the configurations become small with 
increasing the maximum density. 

3.2.2 Relations between physical quantities and maximum 
density 

The relation between the physical quantities and the maxi- 
mum density is an important information and it sometimes 
helps us to understand the stability of equilibrium configura- 
tion. Figure [T2] shows maximum density and baryon mass re- 
lations in the magnetized rotating stars for the parameters, 
a = 12 and p, = 0.1 for the different EOSs. Long-dashed, 
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short-dashed, and dotted Hnes represent the relations for 
q = 0.90, q — 0.80, and q — 0.70, respectively. As a refer- 
ence relation, we also show one in the spherical configuration 
as g = 1.0. All sequences reach the mass-shedding limit at 
q ~ 0.60. Closed circles correspond to the maximum mass 
points. As the axis ratios become small, the maximum mass 
points shift to left, i.e. lower density region. We find that 
the increase of the maximum mass is up to about twenty 
percents for the every EOS employed here. In Figure 1131 
the magnetic fields at the pole are shown as a function of 
the maximum densities in the same sequences. For the se- 
quences with SLy, FPS, and LS EOSs, the magnetic field at 
pole is an increasing function of the maximum density for 
any axis ratio. On the other hand, it is found that only for 
the sequences with Shen EOS, these exist maximum of the 
polar magnetic fields around p ~ 2 x lO^'^g cm^^. 

How about relations between other physical quantities 
and the maximum mass ? Figures 1141 1151 and [TS] are the 
relations between the rotation periods, equatorial radii, and 
maximum magnetic field strength with the maximum den- 
sity. Figure [T^ shows the rotation periods are order of mil- 
lisecond and they decrease with decreasing a value of q. This 
is because these sequences are rotation-dominated type (see 



Table |4| and more centrifugal force is needed to deform the 
equilibrium star largely. These features do not depend on 
the EOSs. From Figure 1151 the equatorial radii are order 
of ten kilometer and decreasing functions of the maximum 
density. We also find equatorial radii increase with decreas- 
ing of q. The maximum magnetic field strength are roughly 
10^'' [G] and are increasing functions of maximum density in 
Figure 1161 Note that their values do not depend on q be- 
cause these sequences belong to rotation-dominated type as 
mentioned above. 

As we have already shown, the magnetic field near the 
rotational axis is like a dipole decaying as Bmax/r^ for r — > 
oo. Now, since we find that Bmax is almost insensitive to 
the values of q with fixing the maximum density. Bp should 
then be proportional to 1/Rp, where Rp is a polar radius. 
For instance in the sequences with SLy EOS, -Rp changes 
only from 11.2[km] to 10.6[km] with decreasing of q from 
0.9 to 0.7. On the other hand, it is found that the change 
in the equatorial radius is much larger up to about twenty 
percents for every EOS employed here when changing q from 
0.9 to 0.7 (see Figure [T^ . For completeness, we finally show 
Af — R relation in Figure 1171 where M and R are baryon 
mass and equatorial radius, respectively. It is noted that the 
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Figure 4. (a) Density contour, (b) magnetic field lines, (c) poloidal component, (d) toroidal component of magnetic field, in the 
meridional plane for the parameters, q = 0.7, a = 16, /t = 0.2, and SLy EOS. The thick, eccentric quarter-ellipse denotes the stellar 
surface. The contours are linearly spaced, i.e., the difference of the physical quantities between two adjacent contours in a tenth of the 
difference between the maximum and minimum values. 



relation is also valid for the spherical stars, implying that a 
hiera rchy of the stiffness becomes so (|Shapiro &i Teukolskvl 
1 19831 ). 

3.3 Model of Proto-Magnetar 

In the previous section, we find that the obtained equilib- 
rium configurations are difficult to be applied for the ob- 
served magnetar. Here we decide to construct an equilib- 
rium configur ation of a proto-magnetar b earing in mind a 
hypothesis bv lDuncan fc ThompsonI (jiggd ) that the magne- 
tars could be originated from the proto-neutron-stars with a 
surface magnetic field of order 10^^ [G] and with its rotation 
period of an order of milliseconds. The models computed 
here should serve as examples of initial models of the hydro- 
dynamic/evolutionary simulations of the proto-magnetar. 

As mentioned in Section 12.21 we can construct a 
model of proto-magnetar with finite temperature by em- 
ploying Shen and LS EOSs. In the actual calculation, the 
maximum density, the parameters associated with mag- 
netic field, axis ratio and entropy per baryon (in unit 
of Boltzmann's constant) are fixed as {pmax,ci, fi,q,s) = 
(lO^'^g cm"^, 20, 0.1, 0.9, 2). Here we take the value of the 



entropy per baryon in the proto-neutron according to 
iPrakash et al.l (|l997h indicating its distributions to be nearly 
uniform as a result of the convection. 

Results are summarized in Table [S] Note that magnetic 
field strength at pole is order of 10^® [G]. Matter and mag- 
netic field structure of these stars are depicted in Figure [18] 
and 1191 As a result, it is found that the features appear- 
ing in the cold magnetized equilibrium configurations above 
are still valid despite of the incursion of the finite temper- 
ature effect on the EOSs. More interestingly, the masses of 
magnetars are also dependent on the EOSs as shown in 
Table [5] They are predicted to be as large as S.OM© for 
the LS EOS, such neutron stars have never been observed 
(ILattimer fc Prakash|[2006l '). Although we have little obser- 
vational information about the masses of magnetar so far, 
our results suggest that we may obtain information about 
the EOS from the observation of the masses of magnetars. 
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Figure 5. Same as FigurcHbut for FPS EOS. 
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Table 4. Physical quantities for the rotating magnetized stars with a = 20, /t = 0.1 and general relativistic 
correction. Note that R is an equatorial radius. 





Bp[G] 


Bc[G] 


mimq] 


P[ms] 


fl[km] 


h 


SLy 


< 


















0.97 


0.663E4 


-17 


0.273E4 


-18 


0.174E401 


0.336E-f01 


0.117E+02 


0.393E+00 


0.90 


0.723E4 


-17 


0.277E4 


-18 


0.185E401 


0.120E+01 


0.124E+02 


0.429E+00 


0.80 


0.826E+17 


0.282E4 


-18 


0.203E401 


0.831E4-00 


0.135E+02 


0.491E+00 


0.70 


0.913E+17 


0.282E4 


-18 


0.222E401 


0.711E4-00 


0.152E+02 


0.560E+00 


0.63 


0.908E4 


-17 


0.267E4 


-18 


0.228E401 


0.684E+00 


0.167E+02 


0.563E+00 


MS 


















FPS 


Q2 < 


















0.97 


0.581E4 


-17 


0.251E4 


-18 


0.135E401 


0.464E4-01 


O.lllE+02 


0.425E4-00 


0.90 


0.635E4 


-17 


0.255E4 


-18 


0.143E401 


0.136E+01 


0.117E+02 


0.459E+00 


0.80 


0.732E4 


-17 


0.259E4 


-18 


0.156E401 


0.924E4-00 


0.128E+02 


0.516E+00 


0.70 


0.813E4 


-17 


0.258E4 


-18 


0.170E401 


0.783E4-00 


0.143E+02 


0.572E+00 


0.63 


0.818E4 


-17 


0.245E4 


-18 


0.176E401 


0.748E4-00 


0.158E+02 


0.574E+00 


MS 


















Shen 


Q2 < 


















0.98 


0.670E4 


-17 


0.277E4 


-18 


0.282E401 


0.372E4-01 


0.140E+02 


0.293E4-00 


0.90 


0.732E4 


-17 


0.280E4 


-18 


0.301E401 


0.106E4-01 


0.149E+02 


0.326E4-00 


0.80 


0.816E4 


-17 


0.282E4 


-18 


0.328E401 


0.761E+00 


0.163E+02 


0.380E+00 


0.70 


0.876E4 


-17 


0.276E4 


-18 


0.353E401 


0.668E+00 


0.183E+02 


0.459E+00 


0.65 


0.855E417 


0.263E4 


-18 


0.358E401 


0.656E4-00 


0.197E+02 


0.463E+00 


MS 


















LS 


Q2 < 


















0.98 


0.613E4 


-17 


0.259E4 


-18 


0.186E401 


0.998E+01 


0.123E+02 


0.333E+00 


0.90 


0.673E4 


-17 


0.263E4 


-18 


0.199E401 


0.125E4-01 


0.131E+02 


0.367E+00 


0.80 


0.762E4 


-17 


0.266E4 


-18 


0.217E401 


0.870E+00 


0.143E+02 


0.419E+00 


0.70 


0.834E4 


-17 


0.263E4 


-18 


0.237E401 


0.747E+00 


0.160E+02 


0.483E+00 


0.62 


0.816E4 


-17 


0.245E4 


-18 


0.244E401 


0.719E+00 


0.180E+02 


0.491E+00 


MS 



















Table 5. Equilibrium configuration of hot neutron star with Shen and LS EOS for {pmax,ci, fi,q, s) = 
(lO^^g cm^**, 20, 0.1, 0.9, 2), where s is the entropy per baryon in unit of Boltzmann's constant. 



q Bp[G] 


Bc[G] 


M[Mq] 


P[ms] 


/?[km] 


h 


Shen 


0.90 6.73E+16 


2.63E417 


1.99E4-01 


1.25E4-00 


13.1E400 


0.367E+00 



LS 



0.90 7.32E4-16 2.80E417 3.01E+01 1.06E+00 14.9E400 0.326E+00 
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Figure 8. (a) Density contour, (b) magnetic field fines, (c) poloidal component of magnetic field, (d) toroidal component of magnetic 
field, in the meridional plane for the parameters, q = 0.7, a = 20, /t = 0.1, and SLy EOS. General relativistic correction is included. The 
thick, eccentric quarter-ellipse denotes the stellar surface. The contours are linearly spaced, i.e., the difference of the physical quantities 
between two adjacent contours in a tenth of the difference between the maximum and minimum values. To make a comparison with the 
Newtonian result easy, these quantities are shown in the non-dimensional unit. 
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Figure 12. Maximum dcnsity-baryon mass relation in the equilibrium configurations with (a) SLy, (b) FPS, (c) Shen, and (d) LS for 
the parameter, a = 20 and (i = 0.1. As a reference, the relation in the spherical configuration are shown as 5 = 1.0. Solid, long-dashed, 
short-dashed, and dotted lines correspond to axis ratio q = 1.0, q = 0.9, q = 0.8, and q = 0.7, respectively. Closed circles shown in the 
each figures are maximum mass points. 
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Figure 13. Magnetic fields strength at pole as functions of the maximum densities in the same sequences shown in Figure [121 
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Figure 14. Rotation periods as functions of tlie maximum densities in tiie same sequences sliown in Figure [12] 
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Figure 15. Equatorial radii as functions of the maximum densities in the same sequences shown in Figure [121 
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Figure 16. Maximum magnetic field strengtli as functions of the maximum densities in the same sequences shown in Figure [T2l 
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Figure 18. (a) Density distribution, (b) magnetic field line (c) poloidal magnetic field (d) toroidal magnetic field in proto-magnetar 
equilibrium configuration with Slion EOS. 
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Figure 19. (a) Density distribution, (b) magnetic field line (c) poloidal magnetic field (d) toroidal magnetic field in proto-magnetar 
equilibrium configuration with LS EOS. 
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4 CONCLUSION 



In this paper, we have investigated equihbrium sequences 
of magnetized rotating stars with fo ur kinds of reaUs- 
tic e quations of state (EO Ss) of SLy (iDouchin fc Haensell 
I2OOII) et al.) FPS (,Pandharipande fc Ravenhalllll989l'). Shen 



( Shen et al.lll998al l. and LS (jLattimer fc Swestvl[l99lh . We 
have empfoyed the Tomimura-Eriguchi scheme to construct 
the equihbrium configurations. 

At first, we have obtained the solution in the regime 
of Newtonian gravity. Basic properties of the magnetized 
Newtonian stars are summarized as foUows: (1) For the 
non-rotating sequences, there exist nearly toroidal config- 
urations, qmin ~ 0, irrespective of the EOSs. The magnetic 
energy stored in the stars increases with the degree of de- 
formation being larger. (2) For the rotating sequences, we 
have categorized the sequences with four kinds of EOSs as 
rotation-dominated (RD) type, magnetic-dominated (MD) 
one, and nearly toroidal one. As a result, the sequences with 
softer EOSs of SLy and FPS are found to belong to RD or 
MD type and the sequences with stiff er EOSs of Shen and 
LS can belong also to the nearly toroidal one. 

(3) We have also focused on the structure of equilibrium 
configurations. Reflecting the stifi'ness of EOSs, the density 
distributions of stars with SLy or FPS EOSs concentrate 
more at the center than those of the stars with Shen or LS 
EOSs. The toroidal fields for SLy and FPS EOSs are found 
to distribute in relatively wider regions in the vicinity of the 
equatorial plane than for Shen and LS EOSs. The poloidal 
magnetic fields are also affected by the toroidal fields be- 
cause poloidal fields are highly distorted where the toroidal 
fields are strong. Regardless of the difference of the EOSs, 
a global configuration of magnetic field line is found to be 
universal, namely the tori of twisted field lines around the 
symmetry axis inside the star and the untwisted poloidal 
field lines, which penetrate the stellar surface to continue to 
the vacuum. 

Then, adding the GR correction to the gravity, we have 
performed the quantitative investigation of the strong mag- 
netized equilibrium stars. As a result, we find that the dif- 
ference due to the EOSs becomes small because all the em- 
ployed EOSs become sufficiently stiff for the large maximum 
density, typically greater than lO^^g cm""^. We have inves- 
tigated the relation between the baryon mass, the magnetic 
field at pole, the rotational period, the equatorial radius, 
and the maximum magnetic field as a function of the maxi- 
mum density. The typical magnetic fields at pole are about 
10^^ G, the periods are about several milliseconds, the radii 
are about ten kilometers and the maximum magnetic fields 
are about lO^^G. The maximum mass is found to be S.QMq 
for SLy EOS, 2.6M0 for FPS EOS, S.SMq for Shen EOS 
and 2.7M0 for LS EOS for g = 0.7 configurations, respec- 
tively. These values are about twenty percents increasing for 
that of the spherical stars. 

Finally, we have computed equilibrium sequences at 
finite temperature for the Shen and LS EOSs aiming 
to construct the equilibrium configurations of the proto- 
magnetars. As a result, it is found that the features ap- 
pearing in the cold magnetized equilibrium configurations 
above are still valid despite of the incursion of the finite 
temperature effect on the EOSs. Since the masses of the 
proto-magnetars are highly dependent on the EOSs, we have 



speculated that one may obtain information about the EOSs 
from the observation of the masses of magnetars. 

It is true that our treatment for the general relativis- 
tic effect is nothing but a crude approximation. Thus we 
consider this study as a prelude to a fully general relativis- 
tic study, however hoping for the moment that our results 
could serve as the initial condition for the hydrodynamic 
evolutionary studies of newly-born magnetars including the 
microphysical EOSs. 
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APPENDIX A: CODE TEST 

To check our code ability, we calcul ate the same sequences 
shown in lYoshida fc Eriguchil (|2006l l with polytropic equa- 
tions of state. We actually obtain two kinds of sequences, 
non-rotating sequence and rotating-sequence. In Table lAll 
some physical quantities of the non-rotating sequence with 
the polytropic index — 3.0, fc = 0.1 and a = 15 (see equa- 



tion 125 



are shown. For the axis ratio g, the 



uppe 



row corresponds to the result by lYoshida fc Eriguchil (|2006l ) 
and the down one do to our result. Table R2l is a result of the 
rotating sequence with polytropic index A'' = 0.5, a = 20, 
and (1 = 0.05. In al l sequences, the rela t ive e rrors for the 
each quantities with lYoshida fc Eriguchil (|2006l Vs result are 
less than one percent and VC is 10"* ~ 10^^. Therefore, we 
confirm that our code works well. 



APPENDIX B: GENERAL RELATIVISTIC 
CORRECTION 

We start from the metric of spherical symmetry space time. 



(Bl) 



Matter is assumed as a perfect fluid, thus the energy mo- 
mentum tensor is 



rpfJ 



P(l + e)4- 



(B2) 



where p, e, and P are a d ensity, specific internal energy , and 
pressure. TOV equation (jOppenhimer fc Volkofj[l939l l IS 

dm . 2 r-, \ 

— 47rr p(l + e), 

2Gm, 



dr 
dP_ 

dr 



(B3) 
(B4) 



where h — 1 + e + P/ p(? is relativistic enthalpy and m is 
mass function defined as 

1 



2Gm 

r 



Equation for the potential ^gr is 
dc&Gfl G ( , 4Trr^P\ 2Gm, 



dr 



G 



(B5) 



This potential, of course, reduces to the gravitational po- 
tential $ in the Newtonian limit. To make clear a general 
relativistic contribution in the potential, we perform the 
Taylor expansion of equation HB5|) with respect to x{r) = 
2Gm{r) /r(? , which is guaranteed to be small than unity in 
star. Result is written in the form. 



dr 



I Anr^pdr' 



-f -!- / Aixr^ pedr 
r ic^r Jq 



(B6) 

n=2 " n=0 

Note that the first term is a contribution from Newtonian 
gravitational potential. So, we define a general relativistic 
correction in the potential as S^gr = c^^gr — $jv. Bound- 
ary condition is derived from a requirement that the outer 



vacuum space time is equivalent to the Schwarzschild space 
time. 



1 / 2GM\ 
$^,.^ln(l-_^ 



(B7) 



where 7i and M are radius and ADM mass of the star. 
From the Taylor expansion with respect to x{R) = 
the boundary condition for the general relativistic correc- 
tion term S^gr is derived. As a result, S^gr is expressed 
in an integral form, 



S^GR{r) 



'(B)" 



n r 

2 °° 



— f 4nf^pedr 





rR ^ °° 

- / df—lATvGf^pe^ \^2;(f) 

OO 



4nGr^P 



(B8) 



This spherically symmetric effective potential is also applied 
in our 2D axisymmetric configuration code, where we first 
compute angular averages of the relevant hydrodynamical 
variables. These are then used to calculate the spherical 
general correction term S^gr{''^) in equation (jBSj. Conse- 
quently, we modify the 2D Newtonian potential (r, 9) to 
obtain the two dimensional general relativistic potential 

*G«(r,e)='i>^''(r,e)+J$j5^(r). (B9) 

The method shown in iMarek et all (|2006l ) has an ambiguity 
in definition of potential and imposes boundary condition 
of their potential at infinity. However, in our method the 
general relativistic correction in the potential is explicitly 
written down and its boundary condition is imposed at the 
stellar surface. So, the two method is a little bit different. 
Verification of our method is shown in the last in this ap- 
pendix. Subsequently, we focus on the Bernoulli equation, 
which is ordinary used to obtain axisymmetric configura- 
tion (e.g. pifl V Under the metric form (|B3|) . four velocity is 
given by 

„M ^ (g-*GH 0,0,0). (BIO) 

The relativistic Bernoulli equation with spherically symmet- 
ric is also written as 



Inh = In It* + C, 



(BU) 



where h and C are relativistic enthalpy, integral constant, 
respectively. Combining equations (|B10|) - HB11|) . we find 



c Inh - 



i^N +5i!>GR) +C. 



(B12) 



Finally we modify the Bernoulli equation of magnetized ro- 
tating case as 



p{u)d 



u + C. 



(B13) 



Note that this treatment may be speculation. To verify our 
treatment, we calculate spherical symmetric stars with four 
kinds of EOS in both ID TOV code and 2D axisymmetric 
code. Relations of mass, both ADM and baryon mass, with 
central density in these stars are displayed in Figure [BT] We 
confirm our treatment works well. 
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Table Al. Physical quantities for the non-rotating sequence with polytropic index N = 3, a = 15, k = 0.1, 
and C l^ = 0. For the axis ratio, the upper row represents Yoshida & Eriguchi's result ijYoshida fc Eriguchil 
|2006|) and the down one does our result. 
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Table A2. Same as Table lAll but for the rotating sequence with polytropic index N = 0.5, a = 20, fc = 0, 1, 
and fl = 0.05. 
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Figure Bl. ADM mass and baryon mass of spherically symmetric stars with SLy, FPS, Shcn, or LS EOSs as a function of central 
density. The solid and dotted lines are the results of TOV code and the cross symbols and filled circles are those of 2D code, respectively. 



